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Abstract. A set of renormalization invariants is constructed using approximate, 
two-flavor, analytic solutions for RGEs. These invariants exhibit explicitly the corre- 
lation between quark flavor mixings and mass ratios from electroweak scale to GUT 
scale in the context of the SM, DHM and MSSM models. The well known empirical 
relations 6*23 oc m^/mfe. On oc md/nii, at electroweak scale can thus be understood 
as the result of renormalization evolution toward the infrared point. The validity 
of this approximation is evaluated by comparing the numerical solutions with the 
analytical approach. It is found that the scale dependence of these quantities for 
general three flavor mixing follows closely these invariants up to the GUT scale. 



I. Introduction 



The mixing of quarks arises from the mismatch between quark mass eigenstates 
and those that participate in the weak interactions. The transformation from the 
weak eigenstates to the mass eigenstates are represented by the CKM matrix, which 
contains three mixing angles and a complex phase, giving rise to the phenomena of 
CP violation. 

A full knowledge of the CKM matrix is therefore crucial in accounting for ob- 
served data and in constructing theories beyond the standard model. Since the 
CKM is obtained from the diagonalization of the mass matrices, it is perhaps not 
surprising that the flavor mixing parameters are correlated to the quark masses. 

It is well known that both the observed quark mass spectrum and the flavor 
mixing parameters exhibit a strong hierarchical structure. The input quark masses 
at Mw scale are roughly taken in their central values as follows [1] (Precisely, as for 
the quark masses, one need consider the QCD effect when the energy is below M^. 
However, in the present context, this is not much relevant for us to draw our main 
conclusion later on) 

rriu ^ 2AMeV,mc ^ 1.27GeV,mt ^ UlGeV 

nid ^ 4.75MeV, m, ^ 0.104GeV, ^ A.2GeV (1.1) 
They satisfy approximately the following relations [2, 3, 4] 

md : nis : mf, oc : : 1 (1-2) 

where A = 0.22. 

There exists a number of parameterizations of the CKM matrix. They include 
the Wolfenstein parameterization [5], as well as the standard parameterization [1], 
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where 



6*13 OC Vub OC A^, 6^23 OC Vcb OC A^. 



(1.4) 



Obviously, both the quark masses and the CKM matrix show the hierarchical struc- 
ture with the parameter A = 0.22, and the CKM mixing angles and the quark mass 
ratios have relations 



at the Mw scale. 

On the other hand, all physical observables remain unchanged under the trans- 
formation of the redefinition of the phases of the quark fields, i.e., the CKM matrix 
can be rephrased as follows 



without changing any physics, where P and Q are diagonal phase matrices. 

In particle physics one of the substantial problems is that of explaining the quark 
masses and their mixings. In the quark sector, both the mass ratios and mixing 
parameters exhibit rather large hierarchies. This observed pattern of fermion masses 
and mixings does not look accidental. In regards to the relations (1.5), when one 
tries to sort out possible clues contained therein, we need also keep in mind that these 
parameters are all measured at low energies. Since the mass matrix evolves with 
energy, it is necessary to bring renormalization effects into the picture. Therefore, 
this relation can be extrapolated to all energies with the aid of renormalization 
group equations. As the energy scale changes, one expects the pattern of regularity 
to evolve as well. 

The problem of the RGE of the quark Yukawa couplings has been studied in 
many papers [6-12]. However, as for the relation Eq.(1.5), one might be tempted to 

ask 

(1) What is the origin of the empirical correlation Eq.(1.5), especially in theory? 

(2) Then what is the correlation between quark flavor mixings and mass ratios 
in all energy scale? 

A complete theory in this regard is certainly lacking and need to be further 
explored. It is the purpose of the paper, without giving a complete rigorous analysis 
of the whole parameter space, to untangle these intriguing questions by using the 
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renormalization equations. We hope that our research work could shed some hght 
on the questions of quark masses and mixings. 

On the other hand, much research work has been focused on the infrared (quasi)fixed 
points of the renormahzation group equations[13-15]. It has been found that, at 
these (quasi)fixed points, the quark masses have a large hierarchical structure, and 
the mixing angles vanish. However, in their neighborhood, small angles will be 
generated through renormalization effects. In this paper, it turns out that, under 
certain simplifying assumptions, one obtains a set of RGE invariants. At low energy 
scale, when these fixed points are approached, due to the existence of these RGE 
invariants, their low energy limits will naturally lead to the correlations described by 
Eq.(1.5). Thus, this suggests that relations between mass ratios and mixing angles 
are dynamical in origin. To assess the accuracy of the simplifying assumptions, we 
compare our results numerically with the full-fledged three flavor RGE. It is found 
that, for a range of parameters, the two methods yield very similar results. 

In this paper we make use of the parameterization of CKM matrix proposed in 
Ref. [16,17]. There, a set of parameters of the CKM matrix was introduced, which 
is independent of the phase matrices P and Q given by Eq.(1.6). In terms of these 
rephrasing invariant parameters, it turns out that the RGE are also simpler than 
those given in terms of other parameters. They are also amenable to a general 
analysis and enable us to draw our conclusions easily in what follows. 

The organization of the paper is as follows. In section II, we start with the 
RGE of the rephrasing invariant parameters of the CKM matrix. A set of RGE 
invariants is then constructed from the approximate analytical solutions assuming 
simple patterns of the CKM matrix. We find, under the infrared approximation, 
they lead to small mixing angles and large mass ratios. Therefore, it shows that, 
the empirical relations for the mixing angles ^13 and ^23 have an origin from the 
renormalization group invariants. The validity and precision of these correlation 
relations are also evaluated by numerical calculations in the presence of three flavor 
mixing with respect to the standard model(SM), the double Higgs model(DHM) and 
the minimal supersymmetric standard model(MSSM) of electroweak interactions 
from low energy up to the high energy scale. It is found that the scale dependence 
of these quantities for general three flavoring mixing follows closely these invariants. 
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with a correction term of the order of O(A^) up to the GUT scale. In addition, 
in order to evaluate the approximation condition for these RGE invariants, their 
renormalization flows as functions of the input mass ratios are also plotted. Finally, 
in subsection 2.4, we present graphically the RGE flow of these quantities running 
from superhigh energy scales to the weak scale, as illustrated in the framework of 
a maximal predictability model [18]. And section III is devoted to discussions and 
summaries. The appendix is devoted to exhibiting the details of the RGE used in 
this paper. 



II. RGE Invariants and Quark Flavor Mixings 



2.1 Rephrasing Invariant Parameters 



There are many possible ways to parameterize the CKM matrix. However, since 
it can be multiplied by rephrasing matrices without changing its physical contents, 
we will choose to use a set which is manifestly rephrasing invariant [16]. This set 
turns out to exhibit the hierarchical feature clearly. In addition, the resulting RGE 
are simpler than those given in terms of other parameters and are amenable to a 
general analysis. These parameters are related to by 
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Also, these {x, y) parameters satisfy the constraints 

{xi + X2 + X3) - {yi + 2/2 + 2/3) = 1 

X1X2 + X2X3 + = yiy2 + 2/22/3 + 2/32/1 (2.1.2) 

which can be derived form the unitarity condition of the CKM matrix with det V = 
1. These constraints leave us with four independent parameters, which is consistent 
with the standard parameterization Eq.(1.3) with three mixing angles and a phase. 
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As a result, all of the measurable quantities \Vijf are directly related to the {x,y) 
variables. These parameters are found to have an explicit hierarchy structure [17] 



x, = 0{l),X2 = 0{X^),xs = 0{X^) 
yi = 0{\'),y, = 0{\').y, = 0{\') 



(2.1.3) 



The evolution equations of the eigenvalues of the Yukawa coupling matrices are 
[6-12] 



,dh. 



lQ7r'^ = hna, + bh^ + 2cY,m^\ 



and 



dt 



(2.1.4) 



(2.1.5) 



Here, t = \n{fi/M]y), where /i is the energy scale and Mw stands for the W bo- 
son mass, i.e. we choose Mw as the renormalization point, and use quark masses 
and mixing angles at to fix the initial conditions. Also, /j and hi denote 
the eigenvalues of the Yukawa coupling matrices for the up-type and down-type 
quarks, respectively. Finally, for the coefficients a„,a^,6, c, we take the notation in 
Ref.[12]. Specifically, the constant c depends on the model used to run the RGE, and 
c = —3/2, 1/2, 1 for SM, DHM, and MSSM respectively. Furthermore, the running 
of the CKM matrix elements has the explicit form 

,dVi^ 



levr^^ = c 
dt 



(2.1.6) 
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In terms of the (x, y) parameters, these equations can be reformulated as follow 
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in which A/^- = — fj, Ahij = hi — h'^, and Ai(A'^ and Bi{B'^ are summarized 
in the appendix. The scahng dependence of the eigenvalue ratios of the Yukawa 
couphngs can be further written as 



lOTT — 

dt 



^ lnr23 ^ 

In rai 
\ lnri2 / 



A/l31 



xi +yi + y2 X2 + y3 



\ 



X2 + 2/2 xi + y-s xs + yi 
\x3 + y3 X2 + yi Xi + y2 ) 



I A/23 \ 



A/: 



31 



V A/12 / 

(2.1.8) 

where = hi/hj, and b = 3,3,6 for SM, DHM, and MSSM respectively. For 
our purposes, it is convenient to reformulate these equations in terms of the sink 
function, 
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The relation Eq.(1.5) can be extrapolated from the weak scale up to the scale Mx by 
means of these renormalization group equations. Using the hierarchy of the Yukawa 
coupling matrices for the up and down type quarks, i.e., ~ 0(1), ~ 0(A^)[9, 
10], as well as the hierarchy of the Xj, yi parameters, Eq.(2.1.7) can be approximated 
by the following equations 
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Likewise, Eq.(2.1.9) can be written as 



dt 



^ ln[sinh(lnr23)] ^ 

ln[sinh(lnr3i)] 
y ln[sinh(lnri2)] j 



/ H2z{xi + yi) H23{x3 + y2) H23{x2 + y^) \ ( A/23 ^ 

^^31 {X2 + 1/2) ^31 {Xl + yz) H31 (X3 + yi) A/31 

\Hu{x3 + y3) Hi2{x2 + yi) Hu{xi+y2) J V A/12 / 



We will now turn to solving these approximate RGE. 



2.2 RGE Invariants 



f2.1.111 



There have been many phenomenological attempts to relate the quark mixing 
parameters to the values of the quark masses. It would be more appealing if the 
correlation is stable against the RG evolution up the high energy scale. Clearly, it 
would be nicer to have the RGE for those observables in analytic form when one 
studies the scale dependence of the CKM elements and the quark mass ratios. In 
fact, as pointed out in Ref.[19], in the neutrino sector, a set of RGE invariants is 
found which bridges the neutrino Yukawa couplings and mixings. There, it is found 
that the infrared fixed point corresponds to the small neutrino mixing angles and 
infinite mass hierarchies. Based on the same footing of quarks and leptons for some 
GUT theories, it motivates us that there might be a plausible scenario in the quark 
sector as well; one would have a natural explanation for large mass hierarchies and 
small CKM mixing angles, if it results from the RGE evolutions as it approaches 
the infrared fixed points. 

We first consider the correlation associated with (2.1.10) and (2.1.11) for the 
simple case of two flavor mixings in (2.1.1). Their analyses are certainly less de- 
manding when the flavor mixing involves the second and third family only, i.e., the 
structure of the CKM matrix has the form 

/lO o\/io o\ 



W 



(2.2.1; 



xi-y2 X2- yi 
V X2-yi xi-y2 J 

and this can be constructed by setting X2 
3^2 ~ 1/3 = in Eq.(2.1.1). On the other hand, in derivation of Eq.(2.1.10) and 



cos 6*23 sin 6*23 
\ sin^ 623 cos^ ^23 / 

0, X3 - ^3 = 0, x3-y2 = 0, and 



1/2 



(2.1.11), we assume there is a very large mass hierarchy between the up and down 
type quarks. As a result, as for the pure mixing between the first two family quarks, 
we do not expect Eqs.(2.1.10, 2.1.11) are good approximations in the first place. 

From Eq.(2.1.10) and (2.1.11), we can find the following renormalization corre- 
lation relations 

-^—tM - + -^—tM -y^) + T, ln[sinh(lnr23)]' = (2.2.2) 
xi — 1/2 dt X2 — yi at at 

or, 

^In[(a;i-y2)(x2-yi)sinh2[ln^])] =0 (2.2.3) 
Hence, we have the RGE invariant quantity 

3^23 = 2[(xi - y2){x2 - yi)Y'^ sinh[ln ^] = const (2.2.4) 

or, in terms of the mixing angle given by Eq. (2.2.1), it can be written as 

^ln(^sin(2e23)sinh[ln^]) = (2.2.5) 

Therefore, for the special structure (2.2.1), the renormalization group invariant has 
the form: 

3^23 = sin(2^23) sinh[ln = const (2.2.6) 

h2 

Obviously, this RGE invariant, which is an approximate analytic solution result- 
ing from Eqs.(2.1.7) and (2.1.9) for two flavor mixing, is invariant for the energy 
evolution from the EW scale up to the GUT scale for the SM, DHM and MSSM re- 
spectively, and the value of the proportionality constant is determined by the physics 
at certain energy scale. 

In the hierarchical limit, i.e., when /13 ^ /i2,we have 

sinh[ln^] f« (2.2.7) 
112 2 /12 

Moreover, at the low energy limit, it is well known that there exists (quasi)fixed 
points for the RGE of the Yukawa coupling matrices [13-15]. At these points, it is 
found that the quarks have zero mixing angles. But in their neighborhood, small 
angles will be generated through renormalization effects. From Eq.(2.1.7), for exam- 
ple, we can easily observe that there is a fixed point at Xi = 1, = j/j = 0, « 7^ 1 for 
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the zero mixing angles. As a result, when these infrared fixed points are approached 
under the low energy approximation, according to Eq. (2.2.6), the small mixing angle 
will naturally lead to the large quark mass hierarchical structure 

sinh[ln^]oc^ (2.2.8) 

Alternatively, it follows that, by assuming the proportionality constant is order 1 at 
certain energy scale and considering Eqs. (2.2.6, 2.2.7), in the neighborhood of these 
(quasi)fixed points we have 

^23 ^ — oc (2.2.9) 

TTlb 

The accuracy of these relations can be evaluated by studying the case of general 
three family mixings in Eq.(2.1.1), i.e., X2 — y2,x^ — |/3,a^3 — Z/2, and X2 — ys, are 
non-vanishing and their contributions in the RG evolution of ^23 are considered. 
3f?23 given by Eq. (2.2.6) becomes scale dependent and we write 

^23 = 2[(xi — y2){x2 — yi)]^''^ sinh[ln ---^] |3x3 = const + correction terms (2.2.10) 

"■2 

Or, approximately, in this three family mixing case, we can get 

sin(2^23) sinh[ln — ] = const + correction terms (2.2.11) 
h2 

We then evaluate 3?23 numerically. In Figure 1, we plot the scale dependence of 
9f?23 for the general three flavor mixings for the SM, DHM, and MSSM, respectively. 
Here, for definiteness, we use the proposed values of the quark masses in (1.1) and 
the initial values of CKM matrix elements in Ref.[l] at t = as an exemplary 
example. Then the RGE invariant 3^23 is given by 3^23 1 t=o = 2[(xi — 1/2) (2^2 — 
yi)]'/'sinh[ln^]|2x2 = 2[(xi -t/2)(x2 -yi)]'/'sinh[ln^]|3x3 = 1-67. Now, 3^23 

ll2 i>'2 

evolves as a function of the momentum from the weak scale all the way up to the 
grand unification scale. For these three models, the numerical correction of the first 
family to the exact analytical solution is explicitly illustrated through the evolution 
trajectories in the graph. 

We find that the variations of 3f?23 are rather slow for each of the models. Specif- 
ically, in the SM, although its variation is relatively fast, the first family interference 
effects are not sizable and 3^23 changes at a level of 0.01% between My/ and the large 
10^^ Gev scale. For the MSSM and the DHM, the derivations from 3f?23 = const are 
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Figure 1: Renormalization evolution of the quantity 3fJ23- The top dotdashed hne 
stands for the SM evolution. The dashed and dotted lines are for the MSSM and 
DHM, respectively. The flat solid line is the RGE invariant for two flavor mixing. 



even less compared to those for the SM. As a result, it is found that the evolution 
of 3^23 closely follows the RGE invariant 3^23 given by the two flavor mixing case. 
The correlation of flavor mixing angles and the mass ratios can be approximated as 
the following 

sin(2^23) sinh[ln ^] = 0(1) + O(A^) (2.2.12) 

with a correction term of the order of O(A^) up to the GUT scale. 

Likewise, consider the mixing of the first and third families only by setting xi — 
y2 = 0,X3 — y2 = 0, X2 — Hi = 0, and X3 — t/i = in Eq.(2.1.1). The structure of the 
CKM matrix has the form 

/ xi - 1/2 X3 - ?/3 \ / cos^6'i::i sin^6'i::i \ 
1 



W 



\ xs-ys Xi-y2 / 



713 u sm t/i3 
1 



(2.2.13) 



\ sin^6'i3 cos'' 6*13 / 



Applying Eq.(2.1.10) and (2.1.11) leads to the following renormalization differenti- 
ation equation 



1 d 

xi — y2 dt 

Namely, 
or 



[xi - 2/2) + ^^^(^3 - ys) + ^ ln[sinh(lnr3i)]' = (2.2.14) 
X3 — ?/3 dt dt 



±ln[{x,-y2){x-s-ys)smh'[ln^])]=0 (2.2.15) 



^ln(^sin(2^i3)sinh[ln^]) = (2.2.16) 



10 



/ 



/ : 

> ^ ' 



5 10 15 20 25 30 



Figure 2: Renormalization evolution of the quantity 3fJi3, using the same notation 
as is Fig.l. 



Thus we have the following RGE invariant corresponding to the exact analytic so- 
lution of Eqs.(2.1.10) and (2.1.11) for the special pattern (2.2.13) 

Ki3 = sin(2^i3) sinh[ln = const (2.2.17) 

All 

Again, the proportionality constant is determined by the physics at certain energy 
scale. And at the low energy limit, we have a small mixing angle near the fixed 
points due to renormalization effect. There, by further assuming the proportionality 
constant is order 1 in Eq.(2.2.17), that gives us 

^^3~^^oc^ocA^ (2.2.18) 

In the general three flavor mixing case, the evolution of 9?i3 becomes energy 
dependent and is evaluated numerically. In Figure 2, we plot the renormalization 
group evolution of 3fti3 for the SM, DHM, and MSSM respectively. 

From Figure 2, one flnds that the difference between the analytic solution and 
the full MSSM is signiflcant compared to other models. But the relative deviations 
for 9f?i3 still less than 0.02% in the whole range of t up to the GUT scale, and the 
DHM is special in that it allows a smallest deviation in the presence of three flavor 
mixings. One can conclude that 

3fJi3 = 2[(xi - 2/2)(x3 - y3)]'/'sinh[ln^] = 0(1) + 0(A6) (2.2.19) 

Ail 

where the constant, as an exemplary example, has been flxed to be 2.88 by consid- 
ering the initial values of the current quark masses and mixing angles at electroweak 
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scale. All in all, 3?i3 is practically constant even up to the GUT scale. The observed 
symmetry (2.2.20) arises as a result of infrared fixed points of the low energy theory. 
In the neighborhood of the fixed point, the mixing angles attain small values and 
the mass hierarchy becomes large. Irrespective of the specific models, the evolution 
of 3?i3 has interesting infrared fixed points. The behavior of the mass and mixing 
hierarchy can thus be explained from the point of view of renormalization group. 



2.3 Mass Hierarchy Dependence 



As described above, a set of renormalization invariants for all energy scales can be 
constructed from Eqs.(2.1.10) and (2.1.11) for special patterns of the CKM matrix. 
They are good approximations of Eqs.(2.1.7) and (2.1.9) under the condition that 
the up and down type quarks have large hierarchical structure. It is interesting 
to investigate the accuracy of these approximations as one varies the degree of 
mass hierarchy. For this purpose we vary the input mass ratios at low energy, and 
calculate, at high energy, the deviations of from constant, as a function of the 
input. For definiteness, we consider such mass dependence for 3^23 only, since the 
behavior of 3fJi3 is similar. 

In Figure 3, we plot the explicit dependence on the up-down type mass ratios 
for the standard model. At the high energy scale, it is found that 

2[(xi - y2)ix2 - yi)Y^^ sinh[ln ^] =const + O(A^) for rrib/mt = 10% 

2[(a;i - y2)ix2 - Vi)]^^^ sinh[ln ^] =const + O(A^) for nib/mt = 20% (2.3.1) 

Explicitly, the approximation of Eq. (2.2.4) is valid up to the order of if the mass 
hierarchy is no more than 10%. And the deviation from Eq. (2.2.4) will be up to A^ 
order when the mass ratio is around 20%. We are thus assured that Eqs.(2.1.10) 
and (2.1.11) are good approximations if the mass ratio is less than 10%. 

Likewise, in Figures 4 and 5, we can find the mass hierarchy dependence of 3?23 
for DHM and MSSM respectively. We use the same notation for the lines as Figure 
3. Thus, for the DHM, 
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Figure 3: Mass hierarchy dependence of the 3?23 renormahzation evolution of the 
SM. The dashed hne is for mt/mi = 20%. The dotted hne is when miy/mt = 10%. 
The dotdashed is the normal SM. The solid line is the RGE invariant for the two 
flavoring mixings. 
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Figure 4: Mass hierarchy dependence of 3^23 renormahzation evolution of DHM. 
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Figure 5: Mass hierarchy dependence of 3fJ23 renormahzation evolution of MSSM. 



SR23 =const + O(A^) for m^lrrit = 10% 
3^23 =const + O(A^) for mb/rrit = 20% 



(2.3.2) 



As for the MSSM, we have 



SR23 = const + O(A^) for rrib/mt = 10% 



= const + 0{X^) for nib/mt = 20% 



(2.3.3) 



Hence, the mass hierarchy dependence of SM and MSSM give us the same approxi- 
mation condition for the vahdity of Eq. (2.2.4), i.e., they both lead to the deviation up 
to order when the mass ratio becomes 20%. However, the DHM is the least sen- 
sitive one, where the deviation is not appreciable even when the mass ratio reaches 
20%. 

Therefore, from the behavior of the evolution of ^23, we can be assured that 
Eq. (2.2.4) is a good approximation for each model if the down-up type mass ratio 
is less than 10%, with a correction term of order A'^. 



2.4 High Energy Scale Evolution 



So far, we have studied the scaling behaviors of physical quantities constructed 
from RGE invariants running from low energy to the grand unification scale. How- 
ever, at the GUT scale, the symmetry group of the model may become larger, and 
this may entail additional symmetries or textures of the quark Yukawa coupling 
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Figure 6: Renormalization evolution of 3?i3 



matrices. Therefore, in order to compare the quark masses and mixings with the 
present experimental data, it is necessary to run the quark masses and mixings from 
the unification scale down to the electroweak scale. Here, it becomes natural and 
instructive to ask the evolution of these physical quantities based on certain sym- 
metry structures from GUT scale to low energy scale. A common phenomenological 
approach is to assume special textures for the quark mass matrices and to derive 
experimentally testable relations among quark masses and mixing angles [18]. We 
focus on the superstring motivated pattern of the quark mass matrices proposed in 
Ref. [18], in which the quark mass matrices have a hierarchy structure with a min- 
imal number of free parameters, constrained by the underlying symmetry theories. 
We then analyze how the relations of Eqs.(2.2.12) and (2.2.21) evolve from 10^^ Gev 
to electroweak scale and examine their behaviors which might impact low energy 
physics. The patterns of the mass matrices have the following explicit hierarchy 
texture 
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After diagonalizing the mass matrixes by unitary transformations, i.e. U^MuU = 
Du,V^MdV = Dd,the CKM matrix is given by Vckm = U^V. Correspondingly, 
the evolutions of quantities 3?i3 and 3^23 are plotted as shown in Fig. 6 and 7. 

For simplicity and definiteness, we have taken the initial input values as a„ = 
120^6?;, ad = 0.9Gev, e = 0.19,w = 1.2 and = 2.2, (f2 = 2.2 + n/2, (^3 = 0. The 
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Figure 7: Renormalization evolution of 3?23- The Solid lines in Figs. 6, 7 are the RGE 
invariants defined by (2.2.6) and (2.2.19) respectively. The dotdashed lines are of 
the SM. The dotted lines are of the DHM, and the dashed lines are for SUSY. 

3?i3 evolution for the SM has the opposite curvature from the other two models, 
which also has the maximum deviation from the central solid line for RGE invariant. 
It reaches around 0.06% in the low energy limit. As for the evolution of 3f?23, all 
the models have the same evolution directions for the RGE fiows, and the maximal 
deviation occurs for the MSSM, at the level of O(A^). 

Our results show that, for mass matrices of the form given in Eq. (2.4.1), renor- 
malization effects are very small, whether we use the SM, DHM or MSSM. If one 
intends to improve the model predictions at low energies, it is necessary to either 
change the input values at the GUT scale, or one could invoke additional symmetry 
groups for the renormalization evolution. Such extra symmetries could arise from 
extra dimension theories [20], for instance. In fact, the extra Kaluza-Klein states in 
these theories could change the RGE trajectories substantially. We hope to return 
to this analysis in the future. 



In this paper we find that the general RGE for the quark mass matrices have 
simple solutions in appropriate two fiavor approximations. These solutions can be 
expressed in the form of energy scale invariants which relate mass ratios and fiavor 



III. Conclusion 
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mixing parameters. Specifically, these invariants are 

2[(xi - y2){x3 - 2/3)]^^^ sinh[ln ^] = const 

til 

2[(xi — y2){x2 — yi)]^^'^ sinh[ln -^j = const (3.1) 

or 

^3 

sin(26'i3) sinh[ln — ] = const 
hi 

sin(2^23) sinh[ln — ] = const (3.2) 
h2 

The validity and precision of these correlation relations are found quite accurate 

when compared with the complete numerical calculations for general three flavor 

mixing from low energy up to the GUT energy scale. It is observed that there is 

little deviation between the two approaches. From the RGE of Eq.(2.1.7) for the 

(x, y) variables, as found before, there is a fixed point at xi = l,Xi = yi = 0,i ^ 1. 

There, physically one can find it will give us zero mixing angles. The physical 

(quasi)fixed points with zero mixing angles for the Yukawa coupling matrices have 

been extensive discussed in the past [13-15]. Close to these points, corresponding to 

nXg nifi 

small mixing angles, Eq.(3.1) implies 6'23 oc — , Oi^ oc — , the well-known empirical 

nib mb 

relations between physical parameters. This suggests a dynamical origin for these 
empirical relations. And the hierarchies in the mass ratios and mixing angles can 
be the result of renormalization evolution as the energy approaches these infrared 
(quasi)fixed points. 

Note that, our approximation depends on the assumed large hierarchy between 
the up and down type quark masses. This works well when applied to mixings 
involving the third family. However, it can not be justified for the mixing between 
the first two light families. Thus, for the Cabibbo angle we do not expect to have 
a relation analogous to Eq.(3.1). Our analysis brings out the qualitative difference 
between the light and heavy families, in agreement with the observed patterns of 
known physical parameters. S 



^ An alternative discussion ean be found in Ref.[21]. There, it is shown that the Cabibbo angle 
and the quark mass ratio have the relation tan \/m77md — constant, which is scale invariant 
from Mz to the Planck scale. 
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We have also analyzed the renormalization evolution starting from high energy 
scales. It is found that, at least for certain existing models, the RGE invariant 
solutions follow closely the numerical calculations. It is our hope that our analysis 
can provide the impetus for further researches along this direction, so that the 
success or failure of a model can be better assessed. The fact that the existence of 
these RGE invariants is rather instructive, especially when one introduces special 
quark texture at high energy scale. It may also put a substantial constraint on 
the parameter space of the quark mass matrices, or may serve as important clues 
in search for new physics symmetry theories. Actually, at high energies, there are 
many possible choices of the physical parameters, in addition to various options 
for the underlying symmetry groups. A more complete analysis of this problem is 
beyond the scope of this work, but will be attempted in a future publication. 
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Appendix 

The A{A') and B{B') in Eq.(2.1.7) are defined as follows: 
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